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Abstract 


The placenta and tumors can exhibit a shared expression profile of proto-oncogenes. The basis 


of placenta-derived heat shock protein gp96, which induces prophylactic and therapeutic T cell 


responses against cancer including hepatocellular carcinoma (HCC), remains unknown. Here, we 


identified the associated polypeptides from human placental gp96 using matrix-assisted laser 


desorption/ionization-time-of-flight and mass spectrometry and analyzed the achieved proteins 


through disease enrichment analysis. We found that placental gp96 binds to numerous peptides 


derived from 73 proteins that could be enriched in multiple cancer types. Epitope-harboring peptides 


from GPC3 and PEGIO were the major antigens mediating anti-HCC T cell immunity. Molecular 


docking analysis showed that the GPC3- and PEG10-derived peptides, mainly obtained from the 


cytotrophoblast layer of the mature placenta, bind to the lumenal channel and client-bound domain 


of the gp96 dimer. Immunization with bone marrow-derived dendritic cells pulsed with recombinant 


gp96-GPC3 or recombinant gp96-PEG10 peptide complex induced specific T cell responses, and T 


cell transfusion led to pronounced growth inhibition of HCC tumors in nude mice. We demonstrated 


that the chaperone gp96 can capture antigenic peptides as an efficient approach for defining tumor 


rejection oncoantigens in the placenta and provide a basis for developing GPC3 and PEG10 


polypeptide-based vaccines against HCC. This study provides insight into the underlying 


mechanism of the antitumor response mediated by embryonic antigens from fetal tissues, and this 


will incite more studies to identify potential tumor rejection antigens from placenta. 


Keywords: Placental gp96, hepatocellular carcinoma, epitopes, GPC3, PEGIO 
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1. Introduction 

Aberrant expression, mutations, and somatic alterations in certain oncogenes and/or tumor 
suppressor genes have been detected in various human malignancies.! The placenta is a 
heterogeneous and invasive tumor-like organ that invades the uterus to enable the exchange of gases, 
nutrients, hormones, and other molecules between the maternal and fetal blood as well as provides 
a barrier that protects the fetus from maternal immune attack.?? Placenta-specific trophoblast cells 
exhibit various similarities with cancer cells such as a rapid proliferation ability, epithelial- 
mesenchymal transition and invasion, cell-cell fusion, induction of angiogenesis, and immune 
escape.^ In addition, the mechanisms regulating the functions of trophoblastic and malignant cells 
are similar. Ferretti et al. ? reported that activation of the phosphatidylinositol phosphoinositide 3- 
kinase/AKT axis is a central feature of signaling pathways to achieve proliferative, migratory, and 
invasive processes in trophoblasts and cancer cells. A recent study reported that human placental 
tissues contain extensive somatic mutations; the mutation rate per trophoblast cluster is similar to 
that in childhood cancers, which are primarily induced by mutations in the uterus. Furthermore, the 
decidual microenvironment that regulates trophoblasts shows many similarities with the tumor 
microenvironment, such as supporting tumor cell survival including immune cells, oxygen levels, 
and energy metabolism." 

The heat shock protein (HSP) gp96, the most abundant chaperone in the endoplasmic reticulum 
(ER), can modulate innate and adaptive antitumor immune responses. As an essential master 
chaperone, gp96 is involved in trafficking of Toll-like receptors and alpha integrins to regulate 
innate immunity as well as T- and B-cell development.*? In terms of antigen-specific immunity, 


tumor-derived gp96-peptide complexes or cell-based gp96—Ig-secreting vaccines elicited specific T 
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cell immunity against parent tumors in both rodent models and clinical trials.!9-!! 


gp96 has been 
proposed to exert multiple roles in T-cell activation, with its unique ability to bind a variety of tumor- 
associated antigens (TAAs) and disease-associated antigens (DAAs) peptides for antigen 
presentation to major histocompatibility complex (MHC) molecules playing a key role in this 
process. ?-? We previously identified a gp96-associated peptide in human liver tumor tissues; this 


16-17 In 


peptide homolog was found to bind to human leukocyte antigen (HLA)-A11 molecules. 
addition, we and others showed that cellular gp96 interacts with MHC-I and transporter associated 
with antigen processing, and gp96 and calreticulin in the ER constitute a relay line for transferring 
associated peptides to MHC-I molecules for T-cell recognition.!*-? These studies demonstrated the 
ability of gp96 to capture and delivery antigens to MHC molecules to generate antigen-specific T- 
cell responses. 

Because of the similarity in antigen expression between cancer and embryonic tissues, gp96 
extracted from the placenta may bind carcinoembryonic antigens and/or proto-oncogene antigens. 
We previously showed that placenta-derived gp96 (Pgp96) can serve as a prophylactic and 
therapeutic vaccine against several cancers, including transplantable melanoma and breast tumors 
in mice.? Moreover, placenta-derived peptides bound to gp96 pulsed dendritic cells (DCs) and 
generated specific CD8+ T cell responses against melanoma and Lewis lung cancer in mouse 
models.?! A clinical trial was recently initiated to test the safety and therapeutic efficiency of 
placental gp96-peptide complexes against non-small cell lung cancer and hepatocellular carcinoma 
(HCC) (No. ChiCTR2100052023). In this study, we isolated gp96-associated peptides from the 


human placenta, identified the T cell epitopes through peptidome analysis using mass spectrometry 


(MS), and explored the mechanism of placental antigen-mediated immunotherapeutic effects 


91 


92 


93 


94 


95 


96 


97 


98 


99 


100 


101 


102 


103 


104 


105 


106 


107 


108 


109 


110 


111 


112 


against HCC. These results may facilitate the design of vaccines and T cell immunotherapies against 


cancer. 


2. Materials and Methods 


Detection and identification of peptides bound to human placental gp96 by mass spectrometry 


Bulk polypeptide fraction was isolated from placental gp96 samples by acid treatment release 


method.'” Prior to mass spectrometry (MS) analysis, polypeptide samples were freeze-dried with a 


centrifuge concentrator, redissolved in 100 uL water, sonicated for 5 min, and centrifuged at 14,000 


Xg to obtain the supernatant. The supernatant was enriched with a polypeptide chip (L121 


mesoporous silicon wafer) and detected using matrix-assisted laser desorption/ionization (MALDI). 


To identify the peptide sequences, the eluents of all samples were combined, lyophilized, and 


redissolved in 15 uL of 0.1% formic acid, and the peptide sequences were identified using a nanoLC- 


Q EXACTIVE (Thermo Fisher Scientific). Analytical separation was performed using gradients of 


HoO/formic acid 100%/0.1% (solvent A) and CH3CN/formic acid 100%/0.1% (solvent B). The 


gradient was run as follows: 0 min 4% B, then to 8% B at 8 min, 22% B at 58 min, 32% B at 70 


min, 90% B at 71 min, 90% B at 78 min at a flow rate of 300 nL/min. The tandem MS scan range 


was 300—1,600 m/z and dynamic exclusion time was 40 s. 


The peptide mixture was analyzed using liquid chromatography-tandem mass spectrometry 


(LC-MS/MS). The SEQUEST HT search engine of Thermo Proteome Discoverer (1.4.0.288) was 


used to search for and identify proteins in the Uniprot-proteome-human (update-20160226) database. 


The search parameters were as follows: no enzyme was required. The initial maximum allowed 


mass deviation of the precursor ion was set to 10 ppm, and the maximum fragment ion mass 
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deviation was set to 20 mDa. Methionine oxidation was used as a variable modification. The filter 


parameter of the results was that the percolator filtered the spectrogram. The delta Cn was «0.1, and 


the false discovery rate was set to 1%. 


Methods 


The mice used, cells culture, the methods of preparation of placental gp96 and recombinant 


gp96 proteins, peptide synthesis and epitope prediction, preparation and immunization strategy of 


pulsed-BMDCs vaccines, animal tumor model and adoptive therapy, enzyme-linked immunosorbent 


spot (ELISPOT) analysis, in vitro tumor cells killing assay, T2-bing assay, tetramer staining and 


flow cytometry, molecular docking analysis, quantitative real-time PCR, Western blot assay and 


statistical analysis are shown in Supplementary Materials and Methods. 


3. Results 


Identification of placental gp96-bound peptides and cytotoxicity testing 


Placenta-derived gp96 (Pgp96) binds intracellular antigenic peptides, including those from 


proto-oncogenes, that may be shared between the placenta and HCC tumors. Therefore, gp96- 


peptide complexes were isolated from the human placenta. The Pgp96-associated peptides were 


extracted, freeze-dried, enriched using L121 mesoporous silicon wafers, and subjected to matrix- 


assisted laser desorption/ionization-time-of-flight (MALDI-TOF) analysis. The sequences of 


peptides were detected using Q-Exactive mass spectrometry (MS). We obtained 1130 and 1107 


peptides with high confidence levels in two individual experiments. The resulting tandem MS data 


of the common peptides were searched in the Uniprot-proteome-human (update-20160226) 
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database using the SEQUEST HT search engine of Thermo Proteome Discoverer (1.4.0.288). We 


identified 73 common proteins in both experiments. Disease enrichment analysis using data from 


the Disease Ontology database demonstrated that these protein genes were enriched in multiple 


cancer types, including liver, breast, and pancreatic cancers and so on (Figure 1A), suggesting that 


the peptide repertoire from placental gp96 covers a large quantity of tumor antigens. 


Then, we examined whether human placental gp96-bound peptides induce tumor-specific T 


cell responses. HLA-A2.1/Kb transgenic mice immunized with Pgp96-pulsed marrow-derived 


dendritic cells (BMDCs) exhibited a significant increase in HepG2-specific T cells, an HLA-A2- 


positive HCC cell line (Figure 1B) and an enhanced cytotoxic effect compared to that in mice 


immunized with recombinant gp96-pulsed or unpulsed BMDCs (Figure 1C). Similar results were 


observed in the human pancreatic adenocarcinoma cell line PANC-1 and lung adenocarcinoma cell 


line A549. In contrast, tumor antigen-free recombinant gp96 expressed in a baculovirus system did 


not induce a tumor specific T cell response, indicating that the antitumor T cell response induced by 


Pgp96 was attributed to the associated tumor antigens. 


We further dissected the HCC antigens from the Pgp96-bound peptides. Through Venn analysis 


using the Gene Expression Omnibus public transcriptome dataset of HCC (GSE101675) and 


placenta (GSE1133), we found that two proteins, GPC3 and PEG10, were present in all three 


datasets (Figure 1D). Another 25 common genes were detected in three other different datasets for 


HCC (GSE101685, GSE105130, and GSE121248), and also included GPC3 and PEGIO (Figure 


1E). 


We analyzed the sources of GPC3 and PEG10 in the placenta. We traced the expression of 


these two genes in single-cell sequencing data of sorted placental cells from first- and second- 
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trimester human placentas.?? The 1567 cells were divided into 13 groups, with each group containing 
specific high-expression genes (Supplementary Figure S1A and S1B), suggesting that the overall 
cell clustering was reliable. Next, we selected the marker genes of different cells based on the results 
reported by Tsang et al. ? and re-annotated the data on the above early placental single cells. As 
shown in Supplemental Figure S1C, re-annotated cells in the 13 groups included nine cell types. By 
searching for cell populations highly expressing GPC3 and PEG10 using the FeaturePlot() function 
of Seurat, GPC3 was found to be highly expressed in stromal and syncytiotrophoblast cells, whereas 
PEG10 was highly expressed in stromal and cytotrophoblast cells (Figure S1D). The GSE173193 
dataset provides 10xGenomics single-cell transcriptome sequencing data for normally delivered 
placentas. As shown in Figure 1F (left panel), the cells were divided into 19 groups after t-stochastic 
neighborhood embedding analysis. Analysis of different placental marker genes showed that 
COLIAI, a marker gene of stromal cells, was not significantly expressed in normal parturient 
placentas, whereas the GPC3 and PEG10 genes were highly expressed in villous cytotrophoblast 
cells annotated by the PARP1 gene (Figure S1B-D, S2, and Figure 1F). The same expression pattern 
was obtained for GPC3 and PEG10 through annotation of the villus cytotrophoblast layer using 
HLA-G- and KRT7- double genes in first- and second-trimester human placentas (Supplementary 
Figure S1B). These results indicate that GPC3- and PEG10-derived peptides bound to gp96 in the 
mature placenta were mainly derived from the cytotrophoblast layer (Figure 1F, right panel). 
Furthermore, we analyzed single-cell RNA sequencing data from HCC tumor specimens based on 
a study by Sun et al. ?* and found that GFP, a marker gene of HCC, was co-expressed with GPC3 
and PEG10 in the same cell (Figure 1G). These results indicate that the expression of GPC3 and 


PEG10 is HCC-specific. 
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Analysis of immunogenicity of placental gp96-bound GPC3 and PEG10 peptides 


Four polypeptides from GPC3 and six polypeptides from PEG10 bound to placental gp96 were 


found. Bioinformatics analysis revealed that these peptides contained abundant HLA-A-restricted 


epitopes (Table 1). Twelve HLA-A*0201-restricted epitopes of PEG10 and eight epitopes of GPC3 


were synthesized. As shown in a T2 binding assay (Figure 2A), three epitopes of GPC3 and four 


epitopes of PEG10 showed high affinity for HLA-A*0201 molecules. Subsequently, HLA-A2.1/Kb 


transgenic mice were subcutaneously vaccinated twice with BMDCs pulsed with a mixture of these 


three peptides of GPC3 or a mixture of the four peptides of PEG10 incorporated with recombinant 


gp96, and peptide-specific CD8* T cells were quantified by analyzing IFN-y* CD8* T cells (Figure 


2B). As shown in Figure 2C, immunized mice showed an activated cytotoxic T lymphocyte (CTL) 


response against epitopes GPC3152-160 and PEG10229.237. The immunogenicity of epitope peptides 


GPC3 152-1600 and PEG10229-237 was confirmed in tetramer (Figure 2D) and enzyme-linked 


immunosorbent spot (ELISPOT) (Figure 2E) assays. Importantly, mice immunized with Pgp96- 


pulsed BMDCs also exhibited specific CTL responses to GPC3152-160 and PEG10220.237, validating 


that placental gp96 was associated with epitope peptides from GPC3 and PEGI0. 


We further modeled the structures of the two defined epitopes from GPC3 and PEGIO bound 


to HLA-A2 using molecular docking analysis. The results showed that both epitopes bind to the 


active site of the a-subunit of HLA-A2 based on the visible hydrogen bonds and strong electrostatic 


interactions. Epitope GPC3152-160 (KVFGNFPKL) interacts with E63, T73, Y99, Q155, Y 159, T163, 


and W167 via hydrogen bonding in the active cavity (Figure 3A). Epitope PEG10229237 


(ALIGQCIHD) interacts with D77, T80, Y84, R97, W147, and Y159 via hydrogen bonding in the 


201 


202 


203 


204 


205 


206 


207 


208 


209 


210 


211 


212 


213 


214 


215 


216 


217 


218 


219 


220 


221 


222 


cavity (Figure 3B). The two epitopes showed low binding energies of -8.0 and -8.8 kcal/mol with 


HLA-A2, indicating highly stable binding. 


We also performed docking modeling of the gp96-dimer with the GPC3 polypeptide containing 


epitope GPC3 152-160 and PEGIO polypeptide containing epitope PEG10229.237 using ZDOCK. As 


shown in Figures 3C and 3D, the red region represents a lumen formed by two loops of the M- 


domain (residues 394—407) and a client protein-binding domain (residues 652-678) of the gp96 


dimer. The bound polypeptides adopted a-helices and/or f-turns as motifs for gp96 dimer 


recognition. GPC3 and PEG10 polypeptides bind near the lumen in different manners. The GPC3 


polypeptide was inserted into the lumen in an a-helical form (Figure 3C), whereas the PEGIO 


polypeptide bound to the pre-Iumenal client protein-binding domain (Figure 3D). As shown in 


Figures 3E and 3F and Table S1, the GPC3 polypeptide can form hydrogen bond interactions with 


the R395, Y678, and D672 residues of gp96, and the PEG10 polypeptide can form hydrogen bond 


interactions with the D672, T675, and S674 residues to achieve stable binding. The binding, 


hydrolysis, and release of ATP can lead to transitions between different conformations of the gp96- 


dimer (Figure S3). Conformational changes in the polypeptide-binding cavity formed by the two 


loops of the M domain and client protein-binding domain, binding state, and binding energy of the 


gp96-dimer with GPC3 or PEGIO polypeptide may also vary (Figure 3G and 3H). Therefore, the 


process of ATP hydrolysis by the gp96-dimer may affect the binding and release of peptides through 


conformational changes. 


GPC3 and PEG10 epitopes within Pgp96-bound peptides induce anti-HCC immune response 


Using the Cancer Cell Line Encyclopedia (CCLE) database, we found that the expression levels 
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of GPC3 and PEG1O were high in HepG2 cells and low in SK-Hep-1 cells (Figure S4). 


Transcriptome analysis, quantitative reverse transcription-polymerase chain reaction (PCR) and 


Western blot assay showed that the expression levels of GPC3 and PEG10 were much higher in 


HepG2 cells than in SK-Hep-1 cells (Figure 4A, 4B and 4C). Twenty-one overlapping genes were 


detected in intersection analysis between protein genes related to placental gp96-bound peptides and 


protein genes that were highly expressed in HepG2 cells compared to in SK-Hep-1 cells, including 


GPC3 and PEG10 (Figure 4D). As shown in Figure 4E, T cells induced by Pgp96, GPC3152-160, Or 


PEG10229.237 epitope exhibited much higher cytotoxicity against GPC3 and PEG10 high-expressing 


HepG2 cells than that against GPC3 and PEG10 low-expressing SK-Hep-1 cells. 


We further investigated the ability of the GPC3 and PEG1O epitopes to inhibit HCC in vivo 


(Figure 4F). Transfer of epitope-induced CD8* T cells resulted in decreased HCC tumor growth and 


weight in BABL/c nu mice compared to the effects of control CD8* T cells (Figure 4G and 4H). 


Similar tumor inhibitory effects were observed in mice treated with placental gp96-induced CD8* 


T cells. 


Polypeptides from GPC3 and PEG10 complexed with recombinant gp96 exhibit anti-HCC 


activity 


As shown in Figure 5A, similar to the above epitopes, these two GPC3 and PEG10 precursor 


peptides complexed with recombinant gp96 exhibited specific cytotoxicity for HepG2 cells but not 


for SK-Hep-1 cells. The ELISPOT assay revealed that GPC3 or PEG10 polypeptide-pulsed BMDC 


immunization resulted in an approximately 6- or 7-fold increase in GPC3 and PEGIO high- 


expressing HepG2 cells-specific T cells compared to control (both P « 0.01). In contrast, no obvious 
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specific T cell response to GPC3 and PEGIO low-expressing SK-Hep-1 cells was detected, 
indicating that the peptide-induced anti-HCC response was GPC3- and PEG10-specific (Figure 5B). 

We further assessed the in vivo therapeutic efficacy of GPC3 and PEGIO polypeptides in 
HepG2-bearing mice. As shown in Figure 5C—D, treatment with CD8* T cells induced by GPC3 or 
PEGI10 polypeptide showed a significant therapeutic effect against HCC compared to that in the 


control (all P « 0.05 or 0.01). 


4. Discussion 

In the current study, using MS-based techniques, we found that abundant tumor-associated 
antigens that were preferentially expressed in multiple tumor tissues bound to placental gp96, and a 
large portion of these antigens were primarily associated with HCC. We further established an MS- 
based immunopeptidomic approach for identifying tumor-specific antigens from the placenta and 
identified GPC3 and PEGIO as HCC antigens with high immunogenicity recognized by T cells 
(Figure 6). These results provide insight into the mechanism of the antitumor response mediated by 
embryonic antigens from fetal tissues or stem cells. 

The structure of gp96 contains four domains: N-terminal, charged linker, middle, and C- 
terminal. Numerous studies showed that exogenous gp96 mediates the uptake and presentation of 
its chaperone-like peptides by antigen-presenting cells to MHC-I molecules and activates specific 
CTLs in vivo.!^ !% 25 However, the structural interaction between gp96 and the associated peptides 
remains unclear. A recent study of the gp96: polypeptide interaction revealed that polypeptides 
associate with the luminal channel that is formed between two loops in the middle domain and 
client-binding domain of the gp96 dimer? In the present study, we showed that placental gp96 
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associates with epitopes-harboring polypeptides derived from GPC3 and PEGIO. The GPC3 and 
PEGIO polypeptides bind to the luminal channel and client-bound domain. Moreover, 
intramolecular conformational changes induced by ATP hydrolysis may regulate the association and 
release of polypeptides from gp96 (Figure 3). Thus, we further uncovered the structural basis of 
gp96 binding of cellular antigenic peptides and validated that the immunopeptidomic analysis of 
gp96-bound peptides are an effective approach for antigen capture and identification. 

An association between placenta and tumor antigens has been widely reported.?7?? In this study, 
we found that placental gp96 was associated with epitopes-harboring polypeptides from GPC3 and 
PEG10 that elicited anti-HCC T cell responses. We further analyzed the source of these antigens in 
the placenta and found that the GPC3 and PEGIO genes were significantly highly expressed in 
villous cytotrophoblast cells of normally delivered placentas (Figure 1F). GPC3 is an oncofetal 
glycoprotein that binds to the cell membrane via a glycosylphosphatidylinositol (GPI) anchor. GPC3 
is expressed in the placenta, numerous embryonic tissues, and various tumors such as HCC but not 
in the healthy adult liver? ovarian carcinoma, and melanomas.???! GPC3 promotes the growth of 
HCC cells in a process involving Wnt/B-catenin signaling and is emerging as a potential therapeutic 
target.? PEG10, a paternally expressed imprinted gene that encodes a cytosolic protein, is primarily 
expressed in the placenta, and its expression levels are elevated in a variety of cancers, including in 
HCC.? PEG10 enhances the proliferation and invasion of cancer cells. In HCC, PEGIO 
overexpression can improve cell invasion ability and decrease cell apoptosis mediated by SIAH 1.?* 
Our current data, along with the above previous studies, clearly demonstrate the similarity of antigen 
expression patterns between the placenta and tumors, and support the design of a new generation of 
proto-oncogene antigen-based cancer vaccines. Further studies are needed to define new placental 
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antigens and determine their potential functions against cancers, including pancreatic cancer, breast 


carcinoma, and lung cancer. 


In summary, we showed that using chaperon gp96 to capture antigenic peptides is an efficient 


approach for identifying tumor rejection antigens generated in the placenta. GPC3 and PEGIO 


antigen-pulsed DCs exhibit therapeutic effects against HCC. Large-scale, systematic studies of 


inherent oncofetal antigens from the human placenta may help reveal novel proto-oncogenes with 


therapeutic potential against cancer. 
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Supplementary figure 1: 
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Analysis and annotation of single-cell transcriptome data of human early 
placentas. Related to Figure 1. (A) Heatmap showing the characteristic genes of different 
cell populations. The t-stochastic neighborhood embedding (t-SNE) algorithm in the R 
package Seurat was used to group cells and draw heat maps of cell specific genes of different 
populations to test the reliability of grouping. (B) Expression of different marker genes. (C) 
Cell population mapping and annotation after 1567 single-cell t-SNE dimensionality 
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442 reduction. (D) GPC3 and PEG10 gene expression in different cell populations in single-cell 


443 transcriptome data. A darker color indicates a higher expression level in the cells. The cell 


444 population circled by the red frame expressed high levels of GPC3 and PEG1O. 
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t-Stochastic neighborhood embedding (t-SNE) dimensionality reduction and gene expression 
of single-cell sequencing data of normal delivery placenta. Related to Figure 1. (A) Expression 
map of different genes visualized using t-SNE dimensionality reduction; red dashed cell group is 
villous cytotrophoblast cells. (B) Expression of different marker genes was visualized using the 
FeaturePlot() function in the Seurat package. A darker color indicates a higher expression level in 


the cells. 


Supplementary figure 3: 


ATP/AMP-PNP bound 
"Twisted V"dimer 


* 90? rotation of N/M 
* N-terminal dimerization 


* $1 swap 
+ ATP 


+ N-term mobile 
domain re-arrangement 


GP96-NMC 
ATP/AMP-PNP/ADP/Apo ; 


ADP bound / ATP/AMP-PNP 
SE dëi -Pi/- ADP bound 
pre -90° rotation of NM £F. Ele aer 
WHAT - N-terminal dimerization 


- $1 swap 


Schematic diagram of the mechanism that ATP hydrolysis of gp96 leads to conformational 


change. Related to Figure 3. In the state of the "twisted V" dimer, gp96's N-terminal domain can 


bind ATP/AMP-PNP (non-hydrolytic), leading to a conformational change in the N-terminal 
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domain and 90? rotation around the N-terminal/middle domain interface, thus promoting N-terminal 
dimerization. At this time, a closed dimer binding ATP/AMP-PNP is formed to prepare for ATP 
hydrolysis. After ATP hydrolysis and ADP release, the transient N-terminal dimer decomposes and 


re-rotates around the N-terminal/middle domain interface to form a mobile N-terminal domain. 


Supplementary figure 4: 
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Heat map of expression levels of GPC3 and PEG10. Related to Figure 4. The expression profile 


data of different hepatocellular carcinoma cell lines were downloaded from the Cancer Cell Line 


Encyclopedia cancer database and analyzed. 


Supplementary Table 1: 


Table S1. Results of molecular docking ZDOCK, discovery studio 2019. 


Receptor Ligand ZDOCK Hydrogen bond Electrostatic Hydrophobic 
score interaction interaction 

gp96- GPCI kos 16.04 A:ARG395:NHI - A:ARG395:NH2 - C:VAL207:CGI - 

dimer C:GLU238:0E2 C:GLU238:0E1 A:TYR652 
A:ARG587:NH2 - A:ARG587:NH1 - C:ARG203 - 
C:ASP204:0D2 C:ASP204:0D1 A:PROS81 
A:GLN668:NE2 - C:VAL219 - 
C:ILE215:0 A:MET662 
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A:TYR678:OH - A:TYR575 - 


C:PHE208:0 C:VAL219 
C:PRO212:CD - A:TRP654 - 
A:TYR678:0H C:ILE215 
gp96- PEGI0O,,,,, 19.6 A:TYR652:0H - A:MET662 - 
dimer C:LEU202:0 C:ILE235 
A:THRÓ7S:N - A:ILE673 - 
C:ILE235:0 C:ILE237 
A:THR675:0G1 - A:TYR652 - 
C:ILE235:0 C:LEU227 
A:TYR677:0H - A:TRP654 - 
C:GLY209:0 C:ILE235 
A:TYR678:0H - A:TYR677 - 
C:ILE231:0 C:LEU210 
C:SER228:0G - A:TYR678 - 
A:GLN578:0 C:CYS234 
A:SER674:CB - 
C:GLN233:0 


Supplementary Materials and Methods: 
Mice 

HLA-A2.1/Kb (HLA-A2) transgenic mice were maintained in the laboratory. BALB/c-nu 
mice were purchased from SPF (Beijing, China). All animals were maintained in a specific 
pathogen-free animal facility and studied at 6-10 weeks of age. Animal experiments were approved 
by the Research Ethics Committee of the Institute of Microbiology, Chinese Academy of Sciences 


(approval number PZIMCAS2011001). 


Cell culture 

HepG2 cells were kindly provided by Prof. X Ye (Institute of Microbiology, Center for 
Biosafety Mega-Science, Chinese Academy of Sciences, Beijing, China). PANC-1 cells were 
cultured in Dulbecco's modified Eagle's medium (HyClone, Logan, UT, USA) supplemented with 


10% fetal bovine serum (FBS). SK-Hep-1, A549 cells and HepG2 cells were cultured in high- 
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glucose Dulbecco's modified Eagle's medium supplemented with 10% FBS. All culture media 


contained 100 U/mL penicillin and 100 ug/mL streptomycin. 


Preparation of placental gp96 and recombinant gp96 proteins 

Placental gp96 proteins (Pgp96) of human were extracted as previously described.”° Briefly, 
after grinding the placental tissue, the supernatant was collected and subjected to ammonium sulfate 
precipitation to obtain the sediment. The dissolved precipitate was purified using ConA-Sepharose 
affinity chromatography, followed by anion exchange chromatography. Endotoxin levels were 
determined using the Limulus Amebocyte Lysate assay («1 EU/mg) (BioWhittaker, Walkersville, 
MD, USA). 

Soluble recombinant human gp96 (Rgp96) was isolated as previously described.** Briefly, 
recombinant human heat shock protein gp96 constructs were subcloned into the pFastBacl vector 
and expressed using the Bac-to-Bac Baculovirus expression system. After filtration and 
concentration, the supernatant was collected, and purified using a HiTrap Q column (GE Healthcare, 
Little Chalfont, UK). After desalination and concentration, the purity of proteins was greater than 


9595, as determined using sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). 


Peptide synthesis and epitope prediction 
The 8- to 10-mer epitope was predicted using long peptides from two websites: SYFPEITHI 


(http://www.syfpeithi.de and ` Immune Epitope Database Analysis ^ Resource 


(http://tools.iedb.org/main/tcell/). Peptides with scores >20 (for SYFPEITHI) or percentile rank 


«0.5 typically have high MHC affinity. All peptides were synthesized by GenScript (Nanjing, China) 
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and were more than 95% pure. 


Preparation and immunization strategy of pulsed-BMDCs vaccines 

Mouse bone marrow-derived DCs (BMDCs) were obtained as previously described.*”-** 
Briefly, bone marrow cells were harvested from the femur and tibia of 5—7-week-old mice, and red 
blood cells were depleted using red blood lysis buffer. The bone marrow cell suspensions were 
cultured in complete RPMI medium containing 10% heat-inactivated FBS, 20 ng/mL granulocyte 
macrophage colony-stimulating factor (GM-CSF), 10 ng/mL interleukin-4 (IL-4), and 50 uM 2- 
mercaptoethanol (Sigma-Aldrich, St. Louis, MO, USA). All cytokines were purchased from 
PeproTech (Rocky Hill, NJ, USA). After 7 days of culture, nonadherent BMDCs were generated 
and collected, and then randomly divided into the following different groups according to the 
experimental design: Pgp96 group: BMDCs incubated with 100 ug/mL placental gp96; Rgp96 group: 
BMDCs incubated with 100 ug/mL recombinant gp96; GPC3 152-150 group: BMDC S treated with 100 
ug/mL recombinant gp96 and 10 ng/mL GPC3 152-160, PEG10229-237 group: BMDC S treated with 100 
ug/mL recombinant gp96 and 10 ug/mL PEG10229237, epitope group: BMDCs treated with 100 
ug/mL recombinant gp96 plus 5 pg/mL GPC3152-160 and 5 pg/mL PEG10229237, positive control (PC) 
group: BMDCs loaded with 50 ug/mL tumor cell lysate; DC group: BMDCs alone. To prepare 
gp96-peptide complexes, gp96 and peptide were mixed and incubated for 10 min at 50?C, followed 
by incubation for 30 min at room temperature. All BMDCs were then cultured in the presence of 1 
ug/mL lipopolysaccharide (LPS) for 24 h at 37?C. After incubation, non-adherent mature BMDCs 
were collected and washed three times with PBS. BMDCs were administered subcutaneously into 
HLA-A2 female mice on days 0 and 7 (5 x 10? cells in 100 uL PBS/mouse). On day 5 of last 
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immunization, the mice were euthanized, and isolated the splenocytes for analyzing by tetramer 
staining, ELISPOT assay and /n vitro tumor killing assay. The same protocol was used in the tumor 


inhibition experiments with polypeptide-induced T cells. 


Animal tumor model and adoptive therapy 

BALB/c nude mice, approximately 4-6 weeks of age, were subcutaneously inoculated with 2x 
10% HepG2 tumor cells in 100 uL of PBS in the right flank. Tumor growth was monitored with 
calipers every other day, and the tumor volume was calculated using the following formula: volume 
= length x width? x 0.5. At four days after tumor cell inoculation, the tumor-bearing mice were 
randomly divided into five groups and immunized intravenously in the tail with 3-4 x 10% T cells 
isolated from the spleens of Pgp96 group mice, Rgp96 group mice, epitope or peptide group mice, 
and PC group mice four times at 3- or 4-day intervals. Adoptive T cells were enriched using a CD8* 
T cell Isolation Kit (Miltenyi Biotec, Gladbach Bergisch, Germany). After treatment, the mice were 


euthanized if the tumor depth exceeded 15 mm. 


Enzyme-linked immunosorbent spot (ELISPOT) analysis 

The IFN-y ELISPOT assays were performed followed by the manufacturer's instructions 
(Mabtech, Mariemont, OH). Briefly, isolated mouse splenocytes were added to the wells of 96-well 
polyvinylidene difluoride plates (BD-Pharmingen, San Diego, CA) in 4-6 replicates and subjected 
to different stimuli. Typical groups included cells with no stimulant, cells with antigens (tumor cell 
lysate or peptides), and cells treated with anti-CD3 (2 ug/mL) and anti-CD28 (1 ug/mL) (PC) and 
incubated at 37?C for the indicated time periods. Spot-forming cells (SFCs) were counted and 
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analyzed using an Immunospot S5 Versa Analyzer (Cellular Technology Limited, Shaker Heights, 


OH, USA). 


In vitro tumor cells killing assay 

Splenocytes were isolated from immunized mice on day 5 after immunization, and then co- 
cultured with BMDCs (at a 1:20 ratio) that had been incubated with the same treatments as 
immunization. Each well also contained 20 IU/mL recombinant murine IL-2 (PeproTech). After 5 
days of culture, splenic effector cells were collected as effector cells. The target cells (HepG2 or 
SK-Hep-1) were labeled with 2.5 uM 5-(6)-carboxy-fluorescein succinimidyl ester (CFSE) dye and 
washed with PBS. We then seeded CFSE-labeled target cells in 100 uL complete RPMI 1640 
medium into 96-well round-bottom plates, and effector cells in 100 uL complete RPMI 1640 
medium were added to the CFSE-labeled target cells at varying ratios for 4—6 h at 37?C. The dead 
cells were labeled with 7-AAD, and cytotoxicity of splenocytes to tumor cells detected using low 


cytometry. 


T2-bing assay 

The T2 cells were incubated with 1 uM human Dim protein (Sigma-Aldrich) and 50 uM 
polypeptide and cultured overnight at 37?C in a 5% CO» incubator. The cells were collected and 
stained with FITC-labeled HLA-A2 antibody, and detected using flow cytometry. The fluorescence 
index (FI) was calculated as described previously;? a fluorescence index >1 was regarded to 


indicate a high-affinity candidate epitope. 
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Tetramer staining and flow cytometry 

A tetramer antibody was prepared according to the instructions of the QuickSwitch™ Quant 
HLA-A‘ 02:01 Tetramer Kit APC (MBL, Nagoya, Japan). The specific operations were as follows. 
First, the lyophilized peptide was dissolved in dimethyl sulfoxide to a final concentration of 10 
mg/mL. Second, 50 uL of QuickSwitch™ Tetramer was mixed with 1 uL peptide in a microtube, 
followed by adding 1 uL Peptide Exchange Factor. These steps were repeated for other peptides. 
Finally, the mixture was incubated for at least 4 h at room temperature in the dark. After the reaction, 
the tetramer antibody was stored at 4?C protected from light until use. Peptide exchange was 
quantified in a Flow Cytometric Sandwich immunoassay according to the manufacturer's 


instructions and analyzed using a Calibur flow cytometer (BD Biosciences, USA) and FlowJo 


Software (Tree Star Inc.). The exchange efficiency of the peptide was 91.9 + 1.8%, indicating that 
the tetramer antibody was prepared successfully. 

Splenocytes obtained from immunized mice were stimulated with 2 ug/mL peptide for 7 days 
at 37°C and 5% CO». The cells were harvested and stained with percp-conjugated anti-mouse CD8, 
FITC-conjugated anti-mouse CD3, and APC-conjugated peptide tetramer antibody. Data were 


acquired using a Fortessa flow cytometer (BD Biosciences, USA). 


Molecular docking analysis 

We used AutodockVina 1.2.2, a computerized protein-ligand docking software, to assess the 
binding energy and interaction patterns between candidate epitopes and HLA-A2 protein.? The 
three-dimensional structures of proteins GPC3 (PDB ID: AF-P51654-F1), PEGIO (PDB ID: 7LGA, 
resolution, 1.9 À), and HLA-A2 (PDB ID: 3HLA, resolution, 2.6 À) were downloaded from the 
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PDB (http://www.rcsb.org/). The molecular structures of the GPC3152-160 and PEG10229.237 epitopes 
were intercepted from them. We first converted all protein and molecular files to PDBQT format, 
removed all water molecules, and added polar hydrogen atoms. The grid box was centered to cover 
the domain of each protein and accommodate free molecular motion. The docking pocket was set 
as a 30 x 30 x 30 Á square pocket, with a grid distance of 0.05 nm for docking. The molecular 
docking results were analyzed, and models were visualized using PyMol 2.5 and LigPlot+ v.2.2. 

AlphaFold2 2.2.4 was used to predict the structure of the GPC3 polypeptide and PEGIO 
polypeptide bound to gp96 identified using MS. The amino acid sequence of human gp96 was 
downloaded from NCBI, and the dimer structure was predicted using the AplhaFold2-multimer 
2.2.4. The ranked O0 structure file was selected as the final simulation structure of the gp96-dimer 
protein by comparing the five result files. ZDOCK version 3.0.2 was used to predict the complex 
model of the receptor and ligand via rigid docking. Interactions between gp96-dimer with GPC3 
polypeptide, or with PEG10 polypeptide were analyzed using Discovery Studio 2019. The final 
results were visualized and analyzed using PyMol 2.5 and LigPlot+ v.2.2. 

We downloaded the gp96 dimer structure in different conformations from the RCSB RBD 
database to analyze its interaction with GPC3 or PEG10 polypeptides under different torsion dimer 
states. The Gp96-NMC/AMP-PNP-bound closed dimer coordinates were obtained from PDB ID 
SULS, GP96-NMC/AMP-PNP-bound "twisted V" dimer coordinates were obtained from PDB ID 
201U, and GP96-NMC/ADP-bound "twisted V" dimer coordinates were obtained from PDB ID 
201V. ZDOCK version 3.0.2 was used to predict the complex model of the receptor and ligand via 


rigid docking. The top 10 ZDOCK scores were analyzed. 
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Quantitative real-time PCR 

Total RNA was extracted from HepG2 or SK-Hep-1 cells using TRIzol reagent (Invitrogen, 
Carlsbad, CA, USA) according to the manufacturer's instructions. RNA was reverse-transcribed 
using PrimeScript™ RT Master Mix to generate cDNA (TaKaRa, Shiga, Japan). cDNA was 
analyzed in real-time PCR using a Rotor-Gene Q (Qiagen, Hilden, Germany). The primers used for 
amplification were as follows: GPC3, (5'-TGGAGAACGTACTGCTTGGTC-3' forward, 5'- 
TCTTCTCAGTTTCAGTGGTGG-3’ reverse); PEG10, (5-GAGAACAGCGGAGAAGGTCC-3' 
forward; 5'-CAAAACCCGCTTATTTCACGC-3' reverse); GAPDH, (5'- 
GGAGCGAGATCCCTCCAAAAT-3' forward; 5'-GGCTGTTGTCATACTTCTCATGG-3' 
reverse). These primers were used with TB Green Premix Ex Taq (TaKaRa). The quantitative PCR 
conditions were used as follows: 95?C for 30 s, followed by 40 cycles of 95?C for 5 s and 60?C for 
45 s. Amplification of specific transcripts was confirmed from the melting curve profiles generated 
at the end of the PCR program. The expression levels of target genes were normalized to that of 


GAPDH and calculated using the comparative cycle threshold (CT) method (24°). 


Western blot assay 
Western blot analysis was performed as previously described.*° Abs against GPC3 (Bioworld 


Technology, catalog no. BS7410), PEG10 (Bioworld Technology, catalog no. MB0036) were used. 


Transcriptome analysis 
Total RNA was prepared from HepG2 or SK-Hep-1 cells dissolved in TRIzol reagent according 
to the manufacturer's instructions. Transcriptome analysis was performed by Novogene Tianjin 
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(Tianjin, China). 


Common database analysis 


All public databases used in this study have been published and can be obtained from the Gene 


Expression Omnibus database of NCBI, including the human tissue-related transcriptome database 


GSE1133; HCC-related transcriptome database GSE101685, GSE105130, and GSE121248; and 


placental single-cell transcriptome database GSE89497. 


Statistical analysis 


All data are expressed as the mean + SD. Statistical variance between groups was analyzed 


using unpaired two-tailed t-test and one-way ANOVA analysis used Graphpad Prism 8 software 


(GraphPad Software, USA). The variance between Kaplan-Meier curves was compared using log- 


rank test. Differences were considered statistically significant at P « 0.05. *P « 0.05, and **P « 


0.01. 


Figure legends 


Figure 1. Immunopeptidome analysis of mass spectrometry-identified peptides from human 


placental gp96. (A) Disease enrichment analysis of proteins from placental gp96-bound peptides 


identified using mass spectrometry. (B-C) Female HLA-A2.1/Kb transgenic mice were immunized 


with BMDCs pulsed with Pgp96, recombinant gp96 (Rgp96), or tumor lysate as a positive control 


(PC), or with unpulsed BMDCs (DC) or PBS alone as a negative control. Splenocytes from 


immunized mice were stimulated with HepG2, PANCI, or A549 whole-cell lysates, and assayed 


using an IFN-y ELISPOT (B). Splenocytes from immunized mice were stimulated and analyzed for 


31 


663 


664 


665 


666 


667 


668 


669 


670 


671 


672 


673 


674 


675 


676 


677 


678 


679 


680 


681 


682 


683 


684 


cytotoxic activity with CFSE-labeled HepG2, PANCI, or A549 cells as target cells (C). All 


experiments were conducted three times. The results are presented as the mean + SD of three 
mice/group. *P < 0.05, **P < 0.01. (D) Venn diagram showing the distribution of shared genes of 
proteins from placental gp96-bound peptides and genes in the Gene Expression Omnibus public 
transcriptome dataset for hepatocellular carcinoma (HCC) (GSE101675) and placenta (GSE1133). 
(E) Venn diagram showing common genes of high-expressed genes in three different datasets for 
HCC. (F-G) t-Stochastic neighborhood embedding dimensionality reduction and gene expression of 


single-cell sequencing data of normal delivery placenta (GSE173193 database) (F) and HCC 2 (G). 


Figure 2. Analysis of HLA-A*0201 restricted epitopes within placental gp96-bound GPC3 or 
PEG10 peptides. (A) Binding affinity of predicted epitopes were quantified using MHC 
stabilization assays with T2 cells using flow cytometry. HBcis.2; and HBcg2-90 peptides served as 
positive (PC) and negative controls (NC). (B) Schema of animal experiments for vaccination 
strategy. Female HLA-A2.1/Kb mice were subcutaneously immunized with BMDCs pulsed with 
the mix of three epitopes of GPC3 or four epitopes of PEG10 using Rgp96 as an adjuvant at weeks 
1 and 2, respectively. (C) Splenocytes were isolated and stimulated with a single epitope (10 ng/mL) 
respectively. Epitopes-specific IFN-y*CD8* T cells were analyzed using flow cytometry. (D-E) 
Female HLA-A2.1/Kb transgenic mice were subcutaneously immunized with BMDCs pulsed with 
Rgp96, Pgp96, Rgp96 complexed with GPC3 152-160 or PEG1022923; epitope, or control BMDCs (DC) 
at weeks 1 and 2. Isolated splenocytes were stimulated and stained with GPC3152-160 or PEG10229. 
237 tetramer. Staining with HBcis-27 tetramer was used as negative tetramer control. The percentage 
of GPC3152-160 or PEG10229-237 epitope-specific tetramer*CD8* T cells were analyzed using flow 
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cytometry (D). Splenocytes were detected in an IFN-y ELISPOT assay. The control peptide 


HBcAgis.2; was used for background evaluation (E). Data are presented as the mean + SD of five 


mice. Ns, not significant; *P < 0.05, **P « 0.01. 


Figure 3. Binding modes of GPC3- and PEG10-derived epitope/HLA-A2 complexes and 


epitope-harboring polypeptide/gp96-dimer complexes. (A-B) Models of GPC3152-160 


epitope/HLA-A2 complexes (A) and PEG10229.237 epitope/HLA-A2 complexes (B) were 


constructed through molecular docking analysis. The three-dimensional structure and interaction 


sites of the binding pocket were visualized using PyMol software (left panel), and two-dimensional 


interactions of epitopes and their targets were exhibited, with hydrogen bonding interactions 


indicated in green (right panel). (C-D) Models of GPC3 polypeptide (aa 149—185)/gp96-dimer 


complexes (C) and PEGIO polypeptide (aa 201—240)/gp96-dimer complexes (D) using ZDOCK. 


The structures of gp96 and GPC3 and PEG10 polypeptide were modeled using AlphaFold2. The left 


panel shows the three-dimensional structure and interaction sites of the binding pocket. The right 


panel represents the binding mode of peptides in binding pocket formed by the two loops of the 


middle-domain and client protein domain of gp96 protein. (E-F) Important binding sites for 


interaction between gp96-dimer and GPC3 (E) or PEG10 (F) polypeptide analyzed using Discovery 


Studio. (G-H) Interaction of gp96-dimer with GPC3 (G) or PEG10 (H) polypeptide under different 


torsion dimer states. The left panel of G or H shows the interaction between GPC3 or PEGIO 


polypeptide with gp96-NMC/AMP-PNP-bound closed dimer, gp96-NMC/AMP-PNP-bound 


"twisted V" dimer, or gp96-NMC/ADP-bound "twisted V" dimer. The right panel of G or H 


represents the ZDOCK score analysis of the top 10 complex conformations. Ns, not significant; *P 
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< 0.05, **P < 0.01, ***P < 0.001, **** P < 0.0001. 


Figure 4. GPC3152-160 or PEG10229-237-specific T cells possess anti-hepatocellular carcinoma 


(HCC) activity. (A) Volcano plots of differentially expressed genes (DEGs) from HepG2 cells 


compared to in SK-Hep-1 cells. DEGs were selected by P < 0.05 and |log 2 (fold-change)| >2. The 


x-axis shows the fold-change in gene expression between HepG2 cells and SK-Hep-1 cells, and the 


y-axis shows the significance of the differences. Colors represent different genes: black indicates 


genes without significantly different expression and red indicates significantly differentially 


expressed genes. (B) Relative expression of GPC3 or PEG10 was analyzed using quantitative real- 


time PCR. Data are the mean + SD of three replicates. (C) Western blot analysis of GPC3 and 


PEG10 expression. (D) Venn diagrams illustrating the number of common genes between placental 


gp96-binding protein genes and genes highly expressed in HepG2 compared to in SK-Hep-1 cells. 


(E) Female HLA-A2.1/Kb transgenic mice were subcutaneously immunized with BMDCs pulsed 


with Rgp96, Pgp96, Rgp96-GPC3 152-160, Rgp96-PEG1 0220-237, or Rgp96-GPC3 152-160 and PEG 10229- 


237 complexes, or unpulsed control DCs at weeks 1 and 2. Five days after immunization, splenocytes 


were harvested for cytotoxic activity with CFSE-labeled HepG2 or SK-Hep-1 cells as target cells. 


(F) Schedule of animal experiments for T cell transfer strategies in HepG2-bearing mice. Mice 


treated with PBS were used as a negative control. (G) Representative images and weights of tumors 


isolated from HepG2-bearing mice. (H) Size of HepG2 tumors was measured at 2-day intervals. n 


= 5 mice/group. PC, positive control. Ns, not significant; *P < 0.05, **P < 0.01. 


Figure 5. GPC3 and PEG10 peptides associated with placental gp96 have therapeutic effects 
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against hepatocellular carcinoma (HCC). Female HLA-A2.1/Kb transgenic mice were 
subcutaneously immunized with BMDCs pulsed with Rgp96-GPC3 polypeptide (aa 149—185), 
Rgp96-PEG10 polypeptide (aa 201-240) complex, or Rgp96 alone as control at weeks 1 and 2. At 
5 days after the last immunization, splenocytes were harvested and stimulated in vitro for 5 days. 
(A) Splenocytes were harvested for cytotoxic activity with CFSE-labeled HepG2 or SK-Hep-1 cells 
as target cells. (B) Splenocytes from immunized mice were stimulated with Hepg2 or SK-Hep-1 
whole cell lysate antigens and assayed in IFN-y ELISPOT assays. (C-D) HepG2 tumor-bearing mice 
(n = 6/group) were intravenously injected with GPC3 or PEG10 polypeptide-induced T cells (3-4 x 
10°/mice) at the indicated times. Representative images and weights of tumors from mice with 


indicated treatment are shown (C). Tumor volumes were measured at 2-day intervals (D). Results 


are the mean + SD. Ns, not significant. Experiments performed twice and showed similar results. 


*P « 0.05, **P « 0.01. 


Figure 6. Schematic illustration of using chaperon gp96 to capture and identify antigenic 


peptides from placenta using mass spectrometry followed by immunopeptidome analysis. 


Placenta-derived gp96 binds abundant tumor-associated antigen peptides associated with multiple 


cancers. Antigenic peptides derived from GPC3 and PEGIO were identified from among gp96- 


bound peptides, which elicit anti-HCC activity. 
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760 Table 1. Prediction of HLA-A-restricted epitopes of placental gp96-bound peptides identified 


761 by mass spectrometry. (Related to Figure 2) 


Gene Peptide sequences (Position) HLA type Epitopes * Position Score" IEDB 
name percentile 
rank * 
GPC3 RDLKVFGNFPKLIMTQVSKSLQV HLA-A*02:01 KVFGNFPKL * 152-160 21 0.09 
TRIFLQALNLGIEV *(149-185) FLQALNLGI 175-183 22 0.21 
RIFLQALNL 173-181 20 12 
HLA-A*11:01 KLIMTQVSK 159-167 21 0.19 
NVLLGLFSTIHDSIQY VOKNAGK HLA-A*02:01 VLLGLFSTI 265-273 26 0.29 
LTTTIGKLCAHSQQRQYRSAY YP TIHDSIQYV 272-280 23 0.05 
EDLFIDKKVLKV (264-321) FIDKKVLKV 313-321 27 0.02 
HLA-A*11:01 STIHDSIQY 271-279 18 0.16 
HLA-A*24:02 AYYPEDLFI 306-314 21 0.06 
HSPLKLLTSMAISVVCFFFLVH HLA-A*02:01 KLLTSMAISV 509-518 25 0.32 
(505-526) 
YTNAMFKNNYPSLTPQAFEFVGE HLA-A*02:01 AMFKNNYPSL 74-83 24 0.77 


FFTDVSL (71-100) 


HLA-A*24:02 NYPSLTPQAF 79-88 23 0.02 
PEGIO AHLATYTEFVPQIPGYQTYPTYA HLA-A*02:01 HLATYTEFV 621-629 25 0.08 
AYPTYPVGFA (620-652) YAAYPTYPV 641-649 21 L 
, HLA-A*11:01 QTYPTYAAY 636-644 11 0.19 
HLA-A*24:02 AYPTYPVGF 643-651 23 0.01 
TYTEFVPQI 624-632 24 0.02 
AAYPTYPVGF 642-651 12 0. 
ALIDQYHEGLSDHIQEELSHLEVA  HLA-A*02:0 ALIDQYHEGL 201-210 26 0.16 
KSLSALIGQCIHIERR! (201-240) ALIGQCIHI 4 229237 24 0.16 
GLSDHIQEEL 209-218 24 0.13 
AM QVQKLTEENTTLREQVEPTPEDE HLA-A*02:0 KLTEENTTL 36-44 24 0.03 
( DDDI (33-59) 
- DHRLVDPHIEMIPGAHSIPSGHVY  HLA-A*02:0 SLSEPEMAAL 548-557 27 0.32 
SLSEPEMAALR (524-558) RLVDPHIEMI 526-535 23 0.32 
NPDMLAPFMAQCQIFMEKSTRD HLA-A*02:0 FMAQCQIFM 96-104 5 0.56 
FSVD (89-114) MLAPFMAQC 92-100 7 0.26 
HLA-A*11:0 AQCQIFMEK 98-106 5 0.15 
HLA-A*24:02 IFMEKSTRDF 102-111 8 0.16 
VRWLSTHDPNITWSTRSIVFDSE HLA-A*02:01 RMYSPIPPSL 474-483 22 0.15 
YCRYHCRMYSPIPPSLPP  (445- WLSTHDPNI 447-455 9 0.55 
485) 
HLA-A*11:01 IVFDSEYCR 462-470 21 0.57 
HLA-A*24:02 MYSPIPPSL 475-483 23 0.01 
RMYSPIPPSL 474-483 11 0.05 
TWSTRSIVF 456-464 12 0.12 


762 " The position of the epitopes according the Protein Group Accessions (P51654-2 for GPC3, ADA087WXK2 for PEG10). The 8- 
763 to 10-mer epitopes were predicted to bind to HLA-A*02 (A*0201), HLA-A*11 (A*1101), and HLA-A*24 (A*2402) by using 


764 different algorithms: the Immune Epitope Database (IEDB) and SYFPEITHI. 


765 : SYFPEITHI (http://www.syfpeithi.de/bin/MHCServer.dll/EpitopePrediction.htm) score of each peptide. A high score indicates 
766 better MHC affinity. 
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112 


° NetMHCpan EL 4.1 (http://tools.iedb.org/main/tcell/) score of each peptide. A low percentile rank indicates better MHC 
affinity. 


s Peptide and epitope were used for anti-HCC animal assay. 
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